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Highly supersaturated electrolyte solutions can be prepared and studied employing an electrodynamic 
levitator trap (ELT) technique. The ELT technique involves containerless suspension of a microdroplet 
thus eliminating dust, dirt, and container walls which normally cause heterogeneous nucleation. This al- 
lows very high supersaturations to be achieved. A theoretical study of the experimental results obtained 
for the water activity in microdroplets of various electrolyte solutions is based on the development of the 
Cahn-Hilliard formalism for electrolyte solutions. In the approach suggested the metastable state for 
electrolyte solutions is described in terms of the conserved order parameter a){j,t) associated with fluc- 
tuations of the mean solute concentration « 0 . Parameters of the corresponding Ginzburg-Landau free 
energy functional which defines the dynamics of metastable state relaxation are determined and ex- 
pressed through the experimentally measured quantities. A correspondence of 96-99 % between theory 
and experiment for all solutions studied was achieved and allowed the determination of an analytical ex- 
pression for the spinodal concentration n 8pin and its calculation for various electrolyte solutions at 298 
K. The assumption that subcritical solute clusters consist of the electrically neutral Bjerrum pairs has 
allowed both analytical and numerical investigation of the number-size N c of nucleation monomers (ag- 
gregates of the Bjerrum pairs) which are elementary units of the solute critical clusters. This has also al- 
lowed estimations for the surface tension a, and equilibrium bulk energy /? per solute molecule in the nu- 
cleation monomers. The dependence of these properties on the temperature T and on the solute concen- 
tration n Q through the entire metastable zone {from saturation concentration n Mt to spinodal n ipin ) is ex- 
amined. It has been demonstrated that there are the following asymptotics: N c = 1 at spinodal concen- 
tration and N c = <x> at saturation. 

PACS number(s): 82.60.Nh, 64.60.My, 82.30.Nr, 82.60.Lf 


I. INTRODUCTION 

Let us consider a binary solution consisting of a strong 
electrolyte and a solvent. Strong electrolytes when dis- 
solved almost completely dissociate into ions. The ther- 
modynamics of these electrolyte solutions are different 
from the thermodynamics of nonelectrolyte binary solu- 
tions. Electrolyte (ionic) solutions possess properties that 
are absent in solutions composed of neutral particles. 
Many of these properties are associated in some way with 
the slow r ~ 1 decay of the Coulomb potential and the cor- 
responding k~ 2 singularity in the Fourier transform. 
The effect of the Coulomb interaction forces between dis- 
solved ions was investigated by Debye-Hiickel [1]. They 
demonstrated that the electrolyte solution energy E ac- 
quires a term proportional to the y power of the mean 
solute concentration n 0 . In binary nonelectrolyte solu- 
tions all thermodynamic functions can be expanded in in- 
teger powers of the concentration n 0 . It is straightfor- 
ward to conclude that the problem of determining ther- 
modynamic quantities of strong electrolyte solutions 
reduces to the problem of a completely ionized gas which 
has been thoroughly discussed in [1] for the case of dilute 
strong electrolyte solutions. 

The physical situation discussed in this paper corre- 
sponds to highly supersaturated strong electrolyte solu- 
tions. Although the problem is of scientific interest and 
practical importance, to the best of our knowledge there 
has been no attempt to provide theoretical insight as well 
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as to develop high-precision measurements. 

It is well known [2-4] that a supersaturated solution of 
a strong electrolyte is in a metastable state. Metastable 
states are characterized by the permanent birth-death 
process of subcritical solute clusters which are homo- 
geneously and randomly distributed over the entire solu- 
tion volume V . For these clusters it is energetically 
favorable to dissolve rather than to grow. This situation 
lasts until the critical or/and supercritical solute clusters 
appear. For these clusters it becomes energetically more 
favorable to grow than to dissolve and nucleation occurs. 
The time required for a nucleation event to occur is 
known as the nucleation induction time. 

Let us consider an electrolyte solution containing n lf 
n 2 y . . *n s ions per unit volume of the various ion species 
denoted by subscripts (1,2, ..., 5 ) which possess the 
charges q\>q 2 , . . .q sy respectively. Using this notation 
the condition of electrical neutrality of the electrolyte 
solution acquires the form 

i iw,-o . (l) 

I = 1 

The solvent in this description is treated as a continuous 
medium, formed from neutral and polar molecules, which 
is characterized only by its dielectric constant e. Since 
the Coulomb forces act between all pairs of ions their 
motion is not entirely random. As a result, the presence 
of an ion at a given location in the solution will affect the 
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spatial distribution of the other ions in its immediate vi- 
cinity. Thus, each ion creates around itself a nonuni- 
formly charged “ion cloud, ” which on average is spheri- 
cally symmetrical. In other words, if we select any par- 
ticular ion, say a j ion, as the center of the solution coor- 
dinate system and consider the density distribution of 
other ions relative to that ion, this density will depend 
only on the distance r from the j ion to the ion con- 
sidered. Thus, let the spatial distribution density n,(r) of 
the i ion in the ion cloud at the spatial point with the ra- 
dius vector r in the 7 -ion-centered coordinate system be 
defined as follows (the subscript j is not attached to the i- 
ion spatial distribution density n,(r) as well as to the elec- 
trostatic potential q> t ( r) at the /-ion location. However, it 
should be kept in mind that these quantities are only 
meaningful in the (moving) /'-ion-centered coordinate sys- 
tem): 

n,= f cf 3 rn,(r) . 

The potential-energy density due to the electronic in- 
teraction between ions and an electric field produced by 
all the others ions around them is 

S 

Uq{t) = — 2 > Qi eZ i 7 ® 

i = l 

where z { is the / ion charge number and e is the unit posi- 
tive charge. In this expression n i 0 is the mean /-ion con- 
centration in solution and <p,( r) is the electrostatic poten- 
tial at the /-ion location produced by all the other ions. 
The electrostatic interaction between ions gives birth, 
therefore, to the following correction E corr (T, V) for the 
electrolyte solution energy E(T,V): 

E(T y V)=E Q (T y V)+E corT (T y V) , 
E m lT,V)=MJ 0 (T,V) , 
where 

A U 0 (T,V)=~^- 2 ’ r = |r| • 

The energy E(T y V) is defined for the case when all ions 
bear their full charges q Xy . . . y q s - The energy E 0 ( T y V) is 
the solution energy when all ions are completely 
discharged and T is the solution temperature. For con- 
venience, we exclude from the E cott (T , V) energy expres- 
sion terms which arise from the ions electrostatic self- 
energy. In these terms the potential (p^ ,• ( r ) is the 
Coulomb field potential of the / ion itself. 

II. THEORETICAL MODEL 
FOR SUPERSATURATED ELECTROLYTE SOLUTIONS 

Supersaturated electrolyte solutions are metastable 
[2-4]. This means that the dissolved solute ions are asso- 
ciated with each other in subcritical solute clusters. 
These clusters appear due to local fluctuations of the 
solute concentration n . It is energetically favorable for 
subcritical solute clusters to dissolve rather than to grow. 
This condition continues until the onset of nucleation. 


when critical and/or supercritical solute clusters form. 
For these clusters it is energetically favorable to grow 
rather than to dissolve. In this way a one-phase metasta- 
ble state of homogeneous supersaturated solution with 
the mean solute concentration n 0 relaxes to a two-phase 
stable state by forming a heterogeneous saturated solu- 
tion with the solute concentration n^. The excess solute 
concentration, An=n 0 -n sat , goes into the creation of 
solid phase. In this two-phase solution the critical and 
supercritical solute clusters are identified with the stable 
solute-rich phase whereas the solution itself is identified 
with the solute-poor phase. The relaxation process from 
the one-phase metastable state to the two-phase stable 
state is known as nucleation. The next stage of the meta- 
stable state relaxation is related to the growth of super- 
critical solute clusters followed by their coalescence. 

In this paper we are interested in the nucleation stage 
of metastable state relaxation in electrolyte solutions. 
The formalism describing this stage of relaxation is based 
on the introduction of the local and time-dependent order 
parameters o^r y t) [3-6]. These order parameters can be 
defined as the scalar functions: 

(u / (r,f) = n I (r, t)—n i0 , 

where rt,(r,f) is the local current concentration of the i 
ions. Therefore, the order parameters w y (r,r) describe 
concentration fluctuations at the spatial-time location 
(r,f). These parameters can be introduced for the both 
phases: the symmetric one, which is associated with the 
undersaturated homogeneous solutions, and the nonsym- 
metric one, which is associated with the saturated- 
supersaturated solutions. These parameters are random 
quantities. Therefore, the experimentally observable 
quantities are their expectation values <u, (r): 

a> t (r)= — f °dta) i (r y t) , 

where t 0 is the macroscopic observation time (in the 
metastable state t 0 is less than the metastable state life- 
time). These expectations are equal to zero: 
6>/ ( r )=0 ( / = 1 , . . . ,s) for the undersaturated states as 
well as for the supersaturated states at the vicinity of 
equilibrium (coexistence or binodal) line. This is under- 
standable since the concentration n ( (r,f) fluctuations 
An^r,/) are negligible in these states and can be of both 
signs with respect to n,(r ,0 with the same probability. 
Therefore, with good accuracy there exist the equalities 
n l (r,/) = n /0 . However, the expectations o),(r) can be 
nonzeros: for the supersaturated states at the 

proximity of spinodal line due to the noticeable and non- 
negligible local fluctuations at the vicinity of the region of 
absolutely unstable states. The theoretical estimation of 
the order of magnitude of the fluctuations An f (r,f) at the 
proximity of spinodal line requires a separate study. 
However, in order to justify the approximation of the 
negligible fluctuations An,(r,r) and, thus, the comparison 
between theory and experiment presented in this paper 
we have taken into account only those metastable states 
in which all the performed measurements had a standard 
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deviation within 5%. 

Since the total solute concentration n as well as the /- 
ion concentrations n i are conserved quantities, the order 
parameters co^Tyt) are conserved order parameters as 
well. This means that the following equalities exist: 


The metastable state in solution is usually created by 
means of a first-order phase transition, for instance by a 
quenching procedure [2-4]. The nucleation stage of the 
metastable state relaxation is described in terms of the or- 
der parameters 6>, (r,U relaxation. It is noteworthy that 
one can consider the order parameter 6),(r,f) relaxation 
independently on the relaxation of other thermodynamic 
quantities only under an assumption that the order pa- 
rameter ^(r,*) relaxation time r y is large enough to 
guarantee the existence of partial equilibrium for every 
specific value of the nonequilibrium order parameter. 
There are two approaches which describe the relaxation 
of the order parameter ^(r,*). They are the Ginzburg- 
Landau formalism developed for the case of noncon- 
served order parameter [3-6,7] and the Cahn-Hilliard 
formalism developed for the case of conserved order pa- 
rameter [3 -6,8,9]. The main assumption in both ap- 
proaches is that nucleation can be associated with a 
Markov’s random process. In our case we will make use 
of the Cahn-Hilliard formalism. Within this formalism 
the nucleation stage of relaxation is given in terms of the 
order parameters relaxation equations: 




2 mw^u) 

a t ' a®,(r,f) ’ ’ 


(4) 


where 

F(a)y<p;t)= f v d*r F(a> t <p;T 9 t) . 

In the above relaxation equations T, is the /-ion mobility 
coefficient and F(co><p;t) is the time-dependent Helmholtz 
free-energy functional for the supersaturated electrolyte 
solution. This functional describes the free energy of the 
order parameter co subsystem and the ion cloud <p subsys- 
tem which interact with each other at the given solution 
temperature and volume. Therefore, the functional 
F(co y <p;t) density F(co,<p;Tyt) should contain the following 
terms: 


F(a)yq>'yTyt)=F((o\Tyt)+F(q)' y Tyt)+F mi (<D,<p\r,t) . (5) 

In expression (5) the Helmholtz free-energy density 
functional F(co;tyt) of the metastable co subsystem can be 
presented as the following Landau expansion with respect 
to even powers of the order parameters a> ( (r,f): 

F(<y;r,t) = 2 

/ = ! 


[Va,(r,f)] 2 ^ a>?(r,r) 


o 


/I, 


1,0 


+ 


Vi 


„ 2 

n i, 0 


( 6 ) 


the correlation radii r c i [5], These radii determine the 
magnitude of distances at which the correlations of the 
order parameters co^tyt) fluctuations decrease 
significantly. The parameters /j, and determine boun- 
daries of the ^-subsystem metastable zone. Since the 
functional F(co\x y t) does not contain interaction terms 
between the various order parameters 
(/ = l,...,s) it is assumed that concentration fluctua- 
tions of the ions of different kinds are independent. 
Therefore, the functional F(co\x y t) represents the simplest 
form for the 6>-subsystem Helmholtz free-energy density 
which allows existence of metastability. 

The ^-subsystem Helmholtz free-energy density func- 
tional F(<p\T y t) has the form 


F(<p;r,r) 


1 

2c 2 


3<p(r,f) 

at 


2 

+ m<p;r,t) 9 


U(<p;T,t) = j[Vq)(Tyt )] 2 , 


(7) 


where c is the light speed in a solution. In this expression 
the term U(<p\x y t) is the ^-subsystem potential-energy 
density. It is assumed that (i) the <p subsystem, which is 
described in terms of the scalar field q>{Xyt) is always in 
equilibrium, and (ii) the ^-subsystem equilibrium state 
cannot be destroyed even by interaction with the co sub- 
system. 

The interaction term F [fiX (cOy<p\Xyt) between the co and 
<p subsystems has the form 

e 5 

F int (<o,<p;r,f) = <p(Tyt)'£ ZjCo^Tyt) , (8) 

E o E f = 1 

where e 0 is the vacuum permittivity constant. Expression 
(8) represents simple electrostatic interaction between the 
ion-cloud potential due to all kinds of ions and an 

electric charge at the spatial-time point (r,/). 

To complete the description of the interacting co and <p 
subsystems one has to supplement Eq. (4), describing the 
evolution of the parameters co,(r,t), with a similar equa- 
tion for the ion-cloud potential <p(r,f). This can be done 
by means of the Euler- Lagrange equation: 


3 &A(co y <p;t) 

3 1 6[3^>(r,f )/3t ] 


I ^ 8A(cOy<p;t) _ S A(cOyq>\t) 
8V<p(r,r) &<p(Tyt) 


(9) 


where A(<Oy<p\t) is the action functional for the supersa- 
turated electrolyte solution. An expression for this func- 
tional in terms of the notations introduced above is 
straightforward: 


A(a)y<p;t)= f d 3 r A(<Oy<p;r y t) , 


A(cOy<p;x,t) = 


2c 2 


3y(r,t) 
3 1 


— U(<P'yTyt) 


~ F(CO\tyt)~ F mX ((Oy<P\Vyt) . 


( 10 ) 


where the parameters B t are proportional to the square of 


One can rewrite Eqs. (4) and (9) as follows: 
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d&iUyt) 

dt 


=r,v 2 


B: . p 

V 2 d>,'(r,r) z,<p( r,f) 

«/, 0 £ 0 £ 


(oA T,t) 

+/L *■’?/ 

n i, 0 


COjilyt) 

3 ' 

Ho 



1 9V( .£ill _ V 2 (r>f )=-g- ^ z<i)j(r,() . 
2c 1 3 1 2 e 0 e , = j 


(11a) 

(lib) 


Let us consider the joint solution of Eqs. (11a) and (lib) 
in the stationary limit. This limit corresponds to the es- 
tablishment of partial equilibrium between solution and 
subcritical ionic clusters. Characteristic times r i 
(/ = 1,2, ... ,s) needed to establish the partial equilibri- 
um should be much less than the induction time. Thus, in 
the stationary limit, i.e., under the conditions that 
da) i (r f t)/dt=0 and dcp(r,t)/dt =0, the spherically 
symmetrical solutions <p(r) and &>,(/*) satisfy the following 
equations: 


B: _ e (Oi(r) 

V 2 6>,-( r) z,<p( r)+ fa 


n 


i, 0 




n 


1,0 


+Vi 


cOjir) 


H 0 


VV(r)= 2 z i co S r ^ * 

e o e i = 1 


= 0 , (12a) 

(12b) 


It is noteworthy that the last equation coincides with the 
electrostatic Poisson equation [1]. 

The properties of the co subsystem which are described 
in terms of the Landau time-dependent Helmholtz free- 
energy density functional F(o>\t,t) are well known near 
the critical point [5,7]. There is a critical point if and 
only if the coefficients /x f - = — | /ij are less than zero 
<0). As it follows from the general order-parameter 
formalism [5,7], the square of the correlation radii r c i 
can be defined in the form r 2 / =^ l -/(2|/xj). Therefore, 
the term (^//iJV^r) is proportional to the ratio 
(r c //, ) 2 , where the length /, characterizes the distance at 
which the field Q)j(r) changes significantly. In order to 
evaluate the ratio r c //, let us assume that the following 
conditions are satisfied. 

(i) The metastable state is far from a critical point, i.e., 
r c « 00 . This is usually the case for electrolyte solutions 
of lower molecular weight solutes. In such solutions the 
critical point has never been observed under normal pres- 
sure. 

(ii) The critical radius of nucleation R c is much greater 
than the correlation radii r ci ( R c »r ci ). The opposite 
case when R c ^r c i corresponds to the unlikely situation 
of immediate nucleation (zero induction time) after the 
solution enters the metastable state. 

Under these conditions the ratio r ci /l t is small: 
r cifh « 1, since R c ~/,. Therefore, Eq. (12a) can be 
simplified as follows: 


<u,(rH 


Vi 

Vi n,;o 


=F,.(r) , 


(13a) 


F i lr)=—^—z i n i0 <p(r) . 
£ o £ M< 


Let us seek for solution of Eq. (13a) in the following form: 


<u,(r) 


2 

m — 0 


Vi 

2 Fi 


m 


(r) . 


(14a) 


The ratios ju, /r/, determine the equilibrium values co i eq of 
the order parameters o),(r): <o i cq /n i0 =((j. i /rf i ) 1/2 , where 
i = 1, . . . , s. The Landau order-parameter formalism for 
the Helmholtz free energy is correct only if the 
Ginzburg-Levanyuk criterion [10,11] is satisfied. This 
criterion requires that the equilibrium values a, eq of the 
order parameters &»,(r) must be much greater than their 
standard deviations: 


— A<u,(r, )A o) i (r 2 )) 

Vi n i,o 

_fij (a> i (r l )6}j(r 2 )) 

Vi [ "?,eq 


From the above criterion it follows that \r) i /2/i i \«l. 
Therefore, considering expansion (14a) it is natural to 
keep only the first two terms which can be found in the 
following form: 


F?(r) 

Vi,o(r)=Fi(r) , { (r )= — — r 

2n i, 0 


(15) 


The next step is to substitute solution for the order pa- 
rameters a)i(r) given by expressions (14a) and (15) into 
the Poisson equation (12b). This substitution leads to the 
following second-order nonlinear ordinary differential 
equation for the scalar field <p{r): 


VV(r)^#e 2 <p(r) + /c 2 <p 3 (r) , 


(13') 


where 


*?=- 


e 

e 0 e 


2 —Ho > 

1 = 1 Fi 


k 2 ~- 


e 

e 0 e 


2 

i=l 


Pi 


Vi*L 0 


It is noteworthy that under the condition k 2 «k \ one can 
reduce the approach developed for supersaturated elec- 
trolyte solutions to the limiting case corresponding to un- 
dersaturated electrolyte solutions. Analytical solution for 
this limiting case was given by Debye and Hiickel [1]. In 
order to obtain the complete correspondence between the 
reduced (k 2 «K|) equation (13') and the Debye-Hiickel 
result for the equation determining the ion-cloud poten- 
tial <p(r) in undersaturated electrolyte solutions one has 
to specify the parameters /i, of the Landau free-energy 
functional (6) as follows: 


fii=fi=-\fi\ , \p\=kT . 


Thus, this specification has allowed us to determine pa- 
rameters /x i of the Landau free-energy functional in the 
temperature region remote from the critical point. 
Therefore, in order to reduce Eq. (13 f ) describing the 
case corresponding to saturated-supersaturated electro- 
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lyte solutions to undersaturated electrolyte solutions it is 
necessary to impose the following constraints on the 
Ginzburg-Landau parameters /i, and 77,: (i) the parame- 
ters Hi are the same for all kinds of ions and are equal to 
“ kTy and (ii) the parameters 17, are all equal to zero 
(17 1 = 0). These constraints are very natural since they 
eliminate the nonlinear terms from the Helmholtz free- 
energy density function F(co)T,t) [see expression (6)] re- 
ducing it to the expression corresponding exactly to the 
Helmholtz free-energy density function of the undersa- 
turated electrolyte solutions [1]. Thus, in terms of Eq. 
03') the effect of the electrolyte solution metastability is 
given by the non-linear term proportional to <p 3 (r). 

Let us seek for solution of Eq. (13') at the /-ion location 
in the following form: 

q>(r)=<pi(r)= £ (c 2 * 2 ) m <p (>I (r) , (14') 

m =0 

where e 2 K 2 is the dimensionless parameter. Substitution 
of this expansion into Eq. (13') allows us to find its first 
terms in the form 


<Pi,o(r)=e 


<p iA (r)=-e 


47re 0 e 
z t A 


A = 


1 + K x a 


E 2 (2k x o) e 
2 K x a 


4 7te 0 e 

/ OO P ^ 

1 d '—’ 

where a is the minimum average distance to which ions, 
both positive and negative, can approach one another. 
Assuming now that the constant e 2 K 2 is much less than 
unit 2 (e 2 /c 2 « 1 ) (This assumption will be verified in Sec. 
VII) one can present the result for the potential q>((r) of 
the supersaturated electrolyte solution as follows: 

12 


q>i(r)=e 


Z: A 


4i re 0 E 


1 —e 2 K~> 


z t A 


47 re 0 e 


' E 2 (lK x a) 
2k x a 


III. THERMODYNAMICS OF THE SUPERSATURATED 
ELECTROLYTE SOLUTIONS 

Let us return to expression (3) for the solution energy 
E ( T, F). It follows from the previous considerations that 
(i) E 0 {T y V) can be identified with the energy of the un- 
dersaturated electrolyte solution consisting of discharged 
ions, and (ii) the energy E(T t V) corrective term 
^corr( F) describes the combined effect due to the elec- 
trostatic interaction between the charged solute ions and 
the solution supersaturation. Thus, substituting into ex- 
pression (3) result (16) for the potential <f>i(r) it is straight- 
forward to obtain the following expression for 
EcoJT t V): 


£ co „(r,K)=Ye 2 K, 


— — e k? 


47re 0 e 
A 


2 z f n i,o 


i = 1 
3 


47TE 0 E 


E 2 (lK x a) ^ 

2 a 2 z i n i.o - 

La i — 1 


(17) 


where the temperature-dependent expressions for k x and 
k 2 are given by relations (13'). 

The Cahn-Hilliard formalism introduced in Sec. II im- 
plies that partial thermodynamic equilibrium exists in a 
metastable state. For the case of nucleation, the lifetime 
of this partial equilibrium can be associated with the in- 
duction time. Therefore, within the induction time all 
equalities of equilibrium thermodynamics exist including 
the relationship between the solution energy corrective 
term E corr ( T f V) and its Helmholtz free energy F(T y V)\ 

E con (T f V) _ 8 F(T,V)-F 0 (T y V) 

t 2 ar t 

In this expression F 0 (T t V) is the Helmholtz free energy 
of undersaturated electrolyte solution consisting of 
discharged ions: 


. (18) 
v 


+o[(e 2 K 2 ) 2 ] . (16) 

It is noteworthy that within the approach developed for 
supersaturated electrolyte solutions the effect of solution 
metastability leads to simple renormalization of the po- 
tential <Pi(r) constant leaving the potential <Pi(r) function- 
al dependence on r the same as it is for the undersaturat- 
ed electrolyte solutions. 


F 0 (T,F) = 6 ) (r,F)+F X [kTrij^n^ + n^V^T, V)] , 

where &( T y V) and 'F,( T, V) (i = 1, . . . ,5) are some con- 
stants at the given temperature and volume. Integration 
of relationship (18) gives the following result for the 
Helmholtz free energy FI T, V): 


F ( r, V) = <?( Ty V ) + V 2 [kTn^n(n^) + ( T y F)/i />0 ] 

1 = 1 


+ 


Ve 2 A 


12t T(kT) W2 (t 0 e) 2 


13/2 


2 z ? n i,o 

i = \ 


Ve*A 2 E 2 (2K x a) 

1 2807r 3 ( /cD 4 ( e 0 e ) 7 a 


s 


£ 

2 ***i,oVi 


2 z Xo 

1 = 1 


i = 1 


(19) 


This expression can be simplified by taking into account 
that the distance a is usually much less than the length 
\/k x (a «\/k x ) which characterizes the dimension of 
the ion cloud of a given ion. Therefore, since qk x «\ 


there exist the approximations 
E 2 (2K x a)^\ , A «1 . 

Using these approximations, the expression for the 1- ion 
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chemical potential derivative dp i {T)/dn it Q in the supersa- 
turated electrolyte solutions acquires the form 


dlij(T) _ i 9 1 F{T y V) 
d"i, o V 9«, 2 o 

2 z . 2 ”/,o 


kT 


z, 4 e 3 


»i,o 16ir(/cD' /2 (e 0 e) 2 


640ir 3 (AcT’) 4 (e 0 e) 7 a 


( 20 ) 


IV. SPINODAL CONCENTRATION FOR THE CASE 
OF BINARY SYMMETRICAL 
ELECTROLYTE SOLUTIONS 

Let us consider the case of binary (i = l,2) symmetrical 
(z, = — z 2 = \z\) electrolyte solutions where the equality 
n i o = n 2 ,o =n o * s true * case Gibbs-Duhem re- 

lationship (which is valid at the states of complete and 
partial equilibriums) acquires the form 

dfi sXv ( Ty Wq ) o dii x (T y n$) d 

N : = — 1 1 1 > 

dn 0 2 an 0 

(21a) 


Analyzing expression (20) for the binary symmetrical un- 
dersaturated electrolyte solutions (t]\— leads to 

the following equality: 

dfi^T) _ dfi 2 (T) _ dfi sh (T,n 0 ) 
dn ,0 dn 20 dn 0 

where |u slt (r,rt 0 ) is the mean solute chemical potential. 
In order to keep this equality for supersaturated electro- 
lyte solutions one has to require that the parameters t) x 
and 7/2 of the ^-subsystem density functional (6) have to 
be the same y]\ — V 2 = V f° r hoth kinds of ions. Under this 
requirement the Gibbs-Duhem relationship for the 
symmetrical supersaturated electrolyte solutions acquires 
the form 


N- 

where 


dfi siv (T,n 0 ) 

dn n 


: —n c 


dMslttr.Kp) 

dn n 


(21b) 


dfi^ t (T,n 0 ) _ kj 

dttn n n 


+ 
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where p sh (T,n 0 ) is the solvent chemical potential and N In this paper we are interested in the solvent activity 
is the number of solvent molecules per unit volume. X %iv (T,nQ) for which there exists the following expression: 


J 


^•siv( 7>o) A slv ( T, 0 )exp 


1 

kTN 



dfxJJ\rn 

dn' 


(22a) 


where A. sIv ( T,0) is the activity of pure solvent. Carrying out integration in expression (22a) we obtain the following rela- 
tion for the activity A. slv ( T,n 0 ): 


k slv (T,n o )=k slv (T,0)cxp 


JL L . 2 1/2 (|z|e) 3 3/2 (zlefrj 2 

N 0 487j-(/cr) 3/2 (E 0 e) 2 ° 128(hr 3 (fcD 5 (e 0 E) 7 a 0 


(22b) 


r 


The contribution due to the electrolyte solution metasta- 
bility comes in the form of the third term in the exponent 
of expression (22b). 

It is well understood [3-6] that at the concentration 
n 0 =n ipi „(T) corresponding to the spinodal line (line 
which separates metastable and unstable states) at the 
given temperature T is where the first derivative 
dfi sit (T,n 0 )/dn 0 of the solute chemical potential 
H iU (T,n 0 ) with respect to the solute concentration n 0 be- 
comes equal to zero: 


dn sh (T,n 0 ) 

dn 0 

Substitution into this equation of result (21b) for 
dfi sh (T,n 0 )/dn 0 allows derivation of the following equa- 
tion for the spinodal concentration n spm ( T): 


= 0 . 


' , 0 = n S pin (r) 


kT+ 


(Uk ) 3 

16ir(2JtD 1/2 (E 0 E) 2 


i wl 
1 spin 


( T) 


( I 2 I e) *V „ . (D=0 

6407r 3 (/cD 4 (e 0 E) 7 a ” sp,n 


(23) 


Solution of this simple square equation with respect to 
n\^(T) and its numerical estimations for various electro- 
lyte solutions will be given in Sec. VII. 


V. SUBCRITICAL SOLUTE CLUSTERS, 

THEIR ASSOCIATION DEGREE, 

AND NUCLEATION MONOMERS 

In the simplest model of electrolyte solutions the dis- 
solved ions can be considered as rigid unpolarized 
spheres contained in a medium of the fixed macroscopic 
dielectric constant e. In this description the nonpolar 
quantum bonds between ions as well as the ion-solvent in- 
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teractions are neglected. Within such a model Bjerrum 
has assumed that two ions of the opposite electric sign at 
the distance r<r ass = (|z|e) 2 /(2e 0 eA:r) are associated 
[12], i.e., from the Bejerrum electrically neutral pair. If 
the ions are of like sign the probability of their associa- 
tion is very small. 

The solution metastable state is characterized by the 
appearance of subcritical solute clusters. Let us assume 
that these clusters are associations of the electrically neu- 
tral Bjerrum pairs. In general, such associations, being 
electrically neutral, may possess dipole moments. How- 
ever, in this paper it is assumed that the dipole moment 
of a subcritical solute cluster, averaged over its lifetime, 
is zero due to the translational and rotational diffusions 
of these clusters. Thus, an interaction between the clus- 
ter dipole moment and the gradient of the ion-cloud po- 
tential, averaged over the cluster lifetime, can be neglect- 
ed. Therefore, within the suggested model the subcritical 
solute clusters can be treated as electrically neutral with 
zero dipole moments associations of the Bjerrum pairs. 
For such clusters it is natural to assume that their associ- 
ated ions do not contribute the correction energy term 
^corr( Tf V) for saturated and supersaturated electrolyte 
solutions [see expression (3)]. Thus, relation (17) for the 
case of saturated and/or supersaturated solutions can be 
presented in the following simple Debye-Hiickel form: 


E mn (T,V)=^-e 

♦ _ e 
K{ t 0 e(kT ) w 2 


2„* 


1 


47re 0 £ 


2 z fK o - 

i = 1 


2 

i = l 


1/2 


(24) 


where n* 0 is the number of i ions which are nonassociated 
into the Bjerrum pairs. 

Let us consider as in the preceding section the case of 
binary symmetrical electrolyte solutions (n^ 0 —n 2i o = n Qy 
7/1 = 7/2 = 17 ). Equating in this case expressions (17) and 
(24) for the solution energy E ( T y V) corrective term 
Eqott^ Ty V) it is straightforward to find n* as the function 
of the solution temperature T and the mean solute con- 
centration n 0 : 


n o ~ n o ( 


= n„ 


n l o /2 V 


l6\ / 2 7r 2 (kT) 7/2 a 


\z\e 

£n£ 


12/3 


(25) 


Therefore, the number n cl (T y n 0 ) of ions associated in the 
Bjerrum pairs can be determined by the following equa- 
tion: 


t _»5<r,» 0 ) , n* 0 (T) 

d(T,n 0 ) = = 1 


= 1 - 


n 'o /2 V 


XeVl^ikT) 1 ^ 


|z|e 

5 

e 0 e 
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(26) 


The ion association degree introduced in this equation, 
d(T y n 0 ) determines the probability that an arbitrary ion 


belongs to the Bjerrum pair. This probability has the lo- 
cal maximum equal to unit [d(T,n 0 )= 1] at the solute 
concentration equal to n 0 max ( T): 


«o,m.*m = 512ir 4 (fc7’) 7 


e 0 e 
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ze 
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(27) 


The temperature-dependent solute concentration 
w o,max(^) corresponds to such a metastable state when all 
ions are associated into the Bjerrum pairs. This means 
that at the concentration «o,max^) the electrolyte solu- 
tion converts into a nonelectrolyte solution since all the 
ions become associated in the electrically neutral Bjerrum 
pairs. As the solute concentration n 0 is increased fur- 
ther, single ions appear until the spinodal concentration 
n g pj n ( T) is reached. 

In the classical approach to statistical mechanics of 
clusters [13] it is assumed that the Helmholtz free energy 
Fi(T f n 0 ) of the cluster containing / particles (in our case 
containing / Bejerrum’s pairs) has the following form: 


F,(T,n 0 )=kT[a(T)l 2/i +pmi-\n(n 0 )l] , (28) 


where n 0 is the mean concentration of solute molecules, 
a( T) is proportional to the surface tension, and /?( T) is 
the equilibrium bulk energy per solute molecule in the 
solute-rich phase (in clusters). Therefore, expression for 
the Helmholtz free energy describes subcritical solute 
clusters as an open (grand) canonical ensemble of the 
Bjerrum pair associations where only the average, rather 
than the exact number of the associated pairs, is specified. 
This means that the solute clusters, the Bjerrum pair as- 
sociations, can exchange the pairs between themselves as 
well as with the solute-poor phase (solution). In such a 
treatment of the Helmholtz free energy the term 
krin(rt 0 ) corresponds to the chemical potential [4,5] of a 
single Bjerrum’s pair in the solution. The grand partition 
function E(T y n 0 ) of the ensemble of subcritical solute 
clusters is equal, by definition, to the ion association de- 
gree d(T y n 0 ) and can be presented as follows: 

E(T,n 0 )=d(T f n 0 )— const X ' 2 ° c ~ F ' iT ' n * )/kT , 

/ = i 


(29) 

where const is the normalization constant and jV c (r,n 0 ) 
is the number size of nucleation monomer which is tem- 
perature and solute concentration dependent. In expres- 
sion (29) we account only for the noninteracting nu- 
cleation monomers (elementary solute units appearing in 
the result of nucleation), aggregations of which constitute 
critical solute clusters. When the number of the Bjerrum 
pairs associated in any subcritical solute cluster reaches 
N c {T y n 0 ) y the nucleation onset begins and nucleation 
monomers form. Equation (29) represent a classical point 
of view for the nucleation of subcritical solute clusters 
[2-6,13]. This equation being combined with result (26) 
for the association degree d(T y n 0 ) y obtained within the 
nonclassical Ginzburg-Landau formalism [3-7], allows 
determination of such classical concepts as the surface 
tension a(T) and the equilibrium bulk energy &(T) as 
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well as the number size N c (T,n 0 ) of nucleation mono- 
mers. Estimation of these concepts will be given in Sec. 
VII. 


VI. EXPERIMENTAL STUDY 
OF SUPERSATURATED ELECTROLYTE SOLUTIONS 

In order to experimentally investigate properties of 
highly supersaturated solutions it is necessary to elimi- 
nate heterogeneous nucleation. For this purpose the elec- 
trodynamic levitation technique has been developed 
[14_16] in order to measure the microdroplet solvent ac- 
tivity A slv (7> 0 ) as a function of the solution temperature 
T and the mean solute concentration n 0 . In this tech- 
nique a micrometer-sized solution droplet is produced by 
the particle generator and after being electrically charged 
it is injected inside the electrodynamic levitator trap 
(ELT), which is incased in a vacuum chamber. The ELT 
consists of three electrodes. In our experiments per- 
formed with the spherical void electromagnetic levitator 
trap (SVELT) [17-21], the top and bottom electrodes are 
spheriods of revolution and the center electrode is torus 
having a spherical cross section. The time-varying volt- 
age V ac (t) is applied to the torus and, relative to the top 
and bottom electrodes interior, a nearly perfect oscillat- 
ing quadrupole potential is produced. This potential, 
even being quite moderate, forces the electrically charged 
solution microdroplet to oscillate vertically along the 
SVELT axis and below its geometric center. Introduc- 
tion of the constant voltage K dc , which is divided equally 
between the top and center, and the center and bottom 
electrodes produce a static interior potential. This poten- 
tial balances gravity at the SVELT null point and, thus, 
allows the microdroplet to be brought to the SVELT 
geometrical center where the ac potential 
K ac U)~c° s (27rW) is zero and the microdroplet comes to 
rest. A diagram of the experimental setup after the solu- 
tion microdroplet being injected is shown schematically 
in Fig. 1. Thus, since the suspended solution microdrop- 
let is containerless (free of foreign surface) the very high 
supersaturations can be achieved before the nucleation 
onset starts. In this way Rubel [14] measured the water 
activity at supersaturations as a function of solute con- 
centration for aqueous phosphoric acid solutions using a 
bihyperboloidal electrodynamic balance trap in a con- 
tinuously flowing system. Using almost the same experi- 
mental setup, Tang and co-workers [15] measured the 
water activity of the supersaturated aqueous electrolyte 
solutions of sodium chloride NaCl and potassium 
chloride KC1. Richardson and Spann [16] used the ELT 
technique to study the water activity of aqueous solutions 
of ammonium sulfate (NH 4 ) 2 S0 4 . In the present study 
we have investigated water activity of the highly supersa- 
turated aqueous electrolyte solutions of sodium bromide 
(NaBr), potassium chloride and bromide (KC1 and KBr), 
and ammonium chloride and sulfate [NH 2 C1 and 
(NH 4 ) 2 S0 4 ]. 

For a stationary microdroplet trapped in the SVELT 
null point its weight mg is balanced by the opposing elec- 
trostatic force qV dc : 



FIG. 1. Schematic drawing of the SVELT experimental ap- 
paratus. 



where q is the microdroplet electrical charge, 2z 0 is the 
distance between the SVELT spheroid electrodes, and C 
is the SVELT geometrical constant. Thus, by means of 
this equation and under an assumption that the micro- 
droplet charge remains unaltered during an experiment 
the relative microdroplet mass changes can be easily 
determined by measuring the balancing dc voltage r d c re- 
quired to balance the weight mg of the charged solution 
microdroplet. The assumption of the constant micro- 
droplet charge q can be easily verified in each experiment. 
This allows us to express the microdroplet solute concen- 
tration n 0 as follows: 


n Q = * — ~ — 1 (in molal units) , 

where (M^) slv is the solvent molecular weight 
[(M^) wate r=18], and and are the balancing 
voltages for the solution (wet) and anhydrous (dry) micro- 
droplets, respectively. 

The two, experimentally justified, assumptions that (i) 
solute is nonvolatile, and (ii) the solution microdroplet is 
in equilibrium with its vapor, allow us to achieve the 
desired mean solute concentration n 0 inside of the solu- 
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tion microdroplet by adjusting the humidity of the 
SVELT chamber vapor. Assumption (ii) means that the 
chemical potentials fi% y (T,n 0 ) and fi^(T,n 0 ) of the sol- 
vent in the liquid (solution) and gaseous (vapor) phases 
are equal. Therefore, the corresponding activities 
^h(T,n 0 )=k^ v (T,0)exp[^ iv (T,n 0 )/kT] and Xg(T, noo ) 
— ktf v (r,0)exp[/j£ v (r,n 0 )/kr] are also equal. Assuming 
that solvent vapor can be considered as an ideal gas it is 
straightforward to demonstrate that [22] 

*4(7>q) _ ^ v (r,n 0 ) _ P(T,n 0 ) 

A4(r, 0) A^r.o) p(t, o) 

_P(T,n 0 ) 

P ( T ) R « (T ’ n ° ) ’ (30) 

where P(T y n 0 ) and R H (T y rt 0 ) are the pressure and the 
relative humidity, respectively, of the solvent vapor 
which is in equilibrium with the solution microdroplet of 
the mean solute concentration n 0 and P( T,0) = / > sat ( T) is 
the pressure of solvent vapor which is in equilibrium with 
the pure liquid solvent [this pressure is the solvent vapor 
saturation pressure P m (T) at the given temperature T]. 
In reality, the solution solvent activity is dependent on 
the microdroplet surface curvature and charge. Howev- 
er, these effects are known to be insignificant for the size 
of the microdroplets being suspended [14]. 

The mean solute concentration n 0 and the correspond- 
ing solution water activity ^ v (7> 0 ) can be easily deter- 
mined by measuring the balancing dc voltage F dc and the 
chamber vapor pressure. Once a microdroplet is caught 
and centered at the SVELT null point, the chamber can 
be evacuated. After the voltage ^dc of the dry micropar- 
ticle is recorded, solvent vapor above the vapor reservoir 
is allowed to bleed back into the chamber until the solid 
microparticle is transformed into solution microdroplet. 
The second evacuation is then commenced at a slower 
rate by adjusting the needle value. This procedure in- 
creases the solute concentration n 0 and thus leads to 
deeper penetration into the supersaturated (metastability) 
region. Evacuation is continued at the slower rate until 
the crystallization point is reached. At that chamber 
pressure the balancing dc voltage K dc drops precipitous- 
ly. After crystallization the evacuation is continued to 
ensure that there has been no charge loss during the cy- 
cle. Therefore, the described experimental procedure al- 
lows continuous recording of the chamber pressure 
P ( T y n 0 ) and the voltage ^dc with the microdroplet at the 
SVELT geometric center. 

As a result of the solvent vapor evaporation, the 
balancing dc voltage V dc decreases steadily and the solu- 
tion microdroplet becomes supersaturated and eventually 
crystallizes. This measurement is repeated several (usual- 
ly three) times to ensure the reproducibility for the given 
solution microdroplet. 

In our experimental study we have investigated the 
SVELT confined microdroplets of aqueous electrolyte 
solutions of sodium bromide (NaBr), potassium chloride 
and bromide (KC1 and KBr), and ammonium chloride 
and sulfate [NH 2 C1 and (NH 4 ) 2 S0 4 ] at the temperature 
T=298 K. The relative humidity R ff (T y n 0 ) has been 


measured versus the mean solute concentration n 0 in the 
wide range An =n 0 — n sat corresponding to the high solu- 
tion supersaturation (riQ^n^y i.e., to the deep solution 
metastable states). Experiments with all solutions were 
performed three times under the same conditions in order 
to guarantee reproducibility. It has been observed that 
the standard deviation cr RH (T,n 0 ) of the relative humidi- 
ty R H (T y n Q ) is a growing function of n 0 . Deep penetra- 
tion into the zone of metastable states was accompanied 
by an increase of a RH (T,n 0 ). This behavior of the stan- 
dard deviation cr RH (T,n 0 ) at high supersaturations is a 
result of local fluctuations of solute concentration and is 
discussed in [23], The increase is conditioned by the cor- 
responding increase of fluctuations of the local solute 
concentration n (r) with approach to the region of unsta- 
ble states (to spinodal line). Therefore, we have retained 
only the R H (T,n 0 ) data where the standard deviation is 
within 5%. 


VII. RESULTS AND CONCLUSIONS 

Figures 2 and 3 show the averaged (over three trials) 
experimental data for the relative vapor humidity 
R H {T y n 0 ) (square points) versus the solute concentration 
n 0 obtained in the SVELT experiment with aqueous mi- 
crodroplets of the ammonium sulfate and ammonium 
chloride. Theoretical treatment (solid line) of these data 
is based on our results presented in Sec. III. According 
to result (21b) for dii s ^(T y n 0 )/dn 0 and relation (30) for 
RH^T y n 0 ) there exists the following expression for 
R h ( T y n 0 ): 


R H (T y n c 


-(M w /\OGO)n 0 WT 7 n Q ) 


where 


\f>(T y n 0 )= 1 + 


(M w ) s] 


48tt( A:T) 3/2 (e 0 e) 2 l 
( \z\e ) 8 7 ? 

mOir\kT) 5 (E 0 e) 7 a 


(M w ) s j v 


In the above expression rfr(T y n 0 ) is the solvent osmotic 
coefficient and the solute concentration n 0 is given in 
molal units [n 0 (in number density units) = (A/^) slv 
Xfl 0 (in molal units)/1000]. The experimental data for 
the relative humidity at the temperature 298 K are extra- 
polated by means of expression (31) in which the solvent 
osmotic coefficient if>(n 0 ) assumes the following form: 


ifj(T y n 0 )— \ + c,( D wq 72 — c 2 ( T)n 0 , 


where 


cAT) = 


2 /2 ( \z \e r 
487r(*7’) 3/2 (£ n e) 2 


c 2 (T)= — 


(M w ),, 


1 280tt 3 ( kT) s ( £ n e ) 1 a 


The numerical expressions for the temperature-dependent 
coefficients c,(D and c 2 (T) can be found by the least- 
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Solute Concentration n „ (in molat units) 


FIG. 2. Dependence of the relative humidi- 
ty R H {T,n 0 ) of the chamber solvent vapor on 
the mean solute concentration ti Q of the 
SVELT confined solution microdroplet. The 
R-H(T>n Q ) measurements correspond to 
7= 298 K. Square points represent experi- 
mental data standard deviation which is within 
5% for each n 0 value. The solid line is the 
theoretical restoration of the obtained experi- 
mental data for supersaturated aqueous solu- 
tion of ammonium sulfate. 


square method applied for the experimentally obtained 
data for Rjj(T,n 0 ) versus n 0 . These constants for each 
system are presented in Table I. The error A (7 = 298 K) 
of the experimental data R H (T t n 0 ) restoration with the 
help of expressions (31) and (32) is very small (within 5%) 
for all the solutions investigated (see Table I). Figures 2 
and 3, for example, give comparison between experiment 
(square points) and theoretical (solid line) restoration by 
means of the results of (31) and (32) for microdroplets of 
supersaturated aqueous solutions of ammonium sulfate 
and ammonium chloride, respectively. In these cases the 
errors A ammon sulf (7=298 K)=4.2929% and 

^amrnon chlor 

(7 = 298 K)= 1.4568%. 

The analytical expressions (32) obtained for the experi- 
mentally measured quantities c^T) and c 2 (7) allow the 
rewriting of Eq. (23) for the spinodal concentration 


^spin 


(7): 

l + ic I (Dn^ 2 n (r)-±c 2 (Df lspin (r)=0. 


Solution of this equation is given by the expression: 

12 


n spi „m=2 + | 


cAT) 


+3 


c 2 (T) 

c,(D 


c 2 m 




cAT) 


cAT) 


1/2 


(33) 


(34) 


Numerical results for the spinodal concentrations 
rt 8pin (7) at 7 = 298 K for various solutions are given in 
Table I. It is noteworthy that in the SVELT performed 
experiments with suspended microdroplets we could 
reach almost spinodal concentration for the supersaturat- 



FIG. 3. Dependence of the relative humidi- 
ty R H {T>n 0 ) of the chamber solvent vapor on 
the mean solute concentration rt 0 of the 
SVELT confined solution microdroplet. The 
Rtf(7,n 0 ) measurements correspond to 
7=298 K. Square points represent experi- 
mental data standard deviation which is within 
5% for each n 0 value. The solid line is the 
theoretical restoration of the obtained experi- 
mental data for supersaturated aqueous solu- 
tion of ammonium chloride. 


Solute Concentration n # (in moial units) 
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TABLE I. (Af^) sk is the solute molecular weight, the concentrations 
molal units. All data correspond to and measurements are obtained at T = 

^sat> ^spim 

298 K. 

and /to. 

max are in 

No. Solute 


z 

^sat 

^spin 

^0,max 

c \ 

c 2 

A (%) 

1 

NaBr 

103 

2 

9.1835 

62.4646 

34.7655 

1.1660 

0.1187 

2.5267 

2 

KC1 

75 

1 

4.7826 

122.5536 

63.9105 

0.5586 

0.0419 

1.7180 

3 

KBr 

116 

1 

5.6962 

66.2428 

34.3240 

0.7355 

0.0753 

0.9711 

4 

NH 2 C1 

54 

1 

7.2733 

34.7944 

15.8927 

0.6152 

0.0926 

1.4568 

5 

(nh 4 ) 2 so 4 

132 

1 

5.7854 

69.6309 

34.2379 

0.5671 

0.0581 

4.2939 


ed aqueous solution of ammonium chloride. For this 
solution the averaged nucleation event occurred at 
n (T = 298 K) = 33.64 (in molal units) whereas the 
analytically derived result for /i spin CT=298 K) = 34.79. 
Therefore, the experimental data obtained for supersa- 
turated aqueous solutions and its treatment within the 
theoretical approach have provided a unique opportunity 
to calculate spinodal concentrations at any temperature 
and, thus, to calculate zones of metastable states. 

In Sec. V we obtained the result that within a region of 
the electrolyte solution metastable states there is a partic- 
ular concentration n 0max (D when the supersaturated 
electrolyte solution experiences a dramatic decrease to 
zero of its electrical conductivity. At this concentration 
the solution becomes nonconductive since all ions become 
associated into electrically neutral Bjerrum’s pairs and 
their associations. Further penetration into metastable 
zone leads back to the appearance of free ions and, thus, 
to nonzero conductivity. Expression (27) for the concen- 
tration n 0 m&x (T) can be rewritten in terms of the experi- 
mentally obtainable coefficients c,( T) and c 2 ( T) y given by 
relations (32), as follows: 


^0,max 


<n= 


3c, (D ] 2 
5 c 2 (T) 


(35) 


Numerical results for the concentrations n 0mas (D at 
7’= 298 K for various solutions are given in Table I. The 


theoretical effect of zero conductivity of the supersatu- 
rate electrolyte solutions at the concentration n 0 nmx (T) 
can be independently verified in the SVELT experiments 
with confined solution microdroplets [24], We are plan- 
ning to perform these experiments in the near future. 
The calculated association degree d(T,rt 0 ) as a function 
of the solute concentration n 0 at T — 298 K is given in 
Fig. 4 for aqueous supersaturated ammonium sulfate 
solution. 

The formalism developed in this paper allows the ex- 
pression of the Landau parameters fi and t/ of the co- 
subsystem free energy in terms of the experimentally 
measured coefficients c , ( T) and c 2 (T)\ 


H=~\kT\ 


iQ7rcu 2 m(E 0 E) 3 (£r) 2 

9c](TK\z\e) 2 


where coefficients c,(D and c 2 (T) are given by expres- 
sions (32). With the help of this result one can justify the 
assumption employed in Sec. II that e 2 K 2 « 1: 

(e 2 K 2 ) mx =A(T)a , 

,, r) _ 2560 ^c 2 ( r)( Eo e) 3 *7>, spin ( T) 

|z|V 

Calculation of the coefficient A( T) results in very small 
values for all systems studied. For example, 

^ammon sulf^ T) » 10 6 . 


1.2 


association degree ( dimensionless ) f d(T ~ 298 K, n 0 ) 



0.2 


FIG. 4. Dependence of the association de- 
gree d{ T,n 0 ) on the mean solute concentration 
n 0 at T = 298 K for supersaturated aqueous 
solution of ammonium sulfate. 
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Phase separation (nucleation) results in the formation 
of the new stable solute-rich phase in the form of critical 
solute nuclei (clusters). In general, these clusters are ag- 
gregates of interacting monomers (nucleation monomers) 
which can be treated as elementary units appeared in the 
result of phase separation. The critical clusters have a 
quite complicated internal structure investigation of 
which requires a separate study. However, estimation of 
the number size of nucleation monomers which we call 
the critical number size N c (T,Kq) can be carried out 
within the approach suggested in this paper. Since the 
direct calculation of the finite sum in expression (29) is 
impossible, it is of use to estimate its lower N cmi „{T,n 0 ) 
and upper N c mMX {T,n 0 ) limits of the critical number size 
N c (T,n 0 ). To perform these estimations we introduce 
the following natural constraints imposed on N c ( T, n 0 ): 

N c (T,n mt (T))=<*> , iV c (7’,n spin (D)= 1 . 

The meaning of these constraints is straightforward: the 
infini tely retarded nucleation at the saturation concentra- 


tion »„,( D requires spontaneous birth of the nucleation 
monomers of infinite number size, whereas at the spino- 
dal concentration /i^l T), where phase separation occurs 
instantaneously, the nucleation monomer number size is 
equal to 1. The last statement means that at the spinodal 
concentration n spin (D the monomers of the new solute- 
rich phase are rather single Bjerrum pairs than their ag- 
gregates. It is well known that at the spinodal concentra- 
tion these monomers (in our case single Bjerrum’s pairs) 
form infinite critical clusters [8,25], The structure of 
these clusters was substantiated to the ramified (fractal) 
[25]. Estimation of jV cmin ( T,n 0 ) implies that the equilib- 
rium bulk energy 0( T)l of nucleation monomer is much 
greater than its surface energy a(D/ 2/3 , i.e., 

f}( T)l l/3 »a(T), whereas estimation of N c mla (T,n 0 ) re- 
quires an assumption that these energies are of the same 
order: fi(T)l x/l ~a{T). Therefore, valuation of relation 
(29) under the above constraints allows us to find expres- 
sions for the lower and upper limits of the critical num- 
ber size N c (T,n 0 ) in the forms 


M c , m JT > n 0 ) = 


1 


N c , dJ T,n 0 ) 


1 


-In 


1- 


ln[£(7> 0 ] 


In 


1 - 


w spin( r > d(T,n 0 ) 
d(T,n spin (T)) n 0 

”s P in( 7 ’) d(T,n 0 ) 


1- 


n Mt (D 


d(T,n spin (T)) n 0 


[1 -$mn 0 ] 


(36a) 


(36b) 


where 

1 l d( r,n spin (D) 

§ n mt (T) n spjn (D d{T,n M (T)) 

In the above expressions n^iT) is the solute concentra- 
tion at saturation (see Table I). The obtained experimen- 
tal results at T— 298 K allow us to investigate how the 
lower N cmin (n 0 )=N c>min (T=298 K, » 0 ) and upper 
N c , m Jn 0 )=N Ctma (T= 298 K, n 0 ) limits of the critical 
number size N c (n 0 ) depend on the solute concentration 
n 0 . This requires substitution into expressions (36a) and 
(36b) of result (34) for n spin (T) and table values for 
n Mt ( T) as well as relation (26) calculated for the given 
n ipt JT) and n M (T). Figure 5 gives the dependence of 

^c,min(«o) and N c, maxK) on «0 in the ran 8 e [ ” spin sat ] 

for the aqueous supersaturated electrolyte solution of am- 
monium sulfate. In Table II there are values for 
JV c Jn 0 ) and N cmxx (n 0 ) at the solute concentrations 


I 

1. 1 and n spin for various solutions. Analysis of Fig. 
5 together with the analogous figures for the other super- 
saturated solutions studied shows that for all the cases 
the critical number size N c {n 0 ) is a decreasing function of 
solute concentration n Q with the following properties: (i) 
N c (n Mt (T))= <x>, (ii) N c (n ip (T))=\, and (iii) N c (T,n 0 ) 
has a local maximum at the concentration fi 0 = n 0 miix ( T). 
Existence of this local maximum can be understood if one 
takes into account that at the concentration T) all 

ions become associated into the electrically neutral Bjer- 
rum pairs and electrolyte solution becomes nonelectro- 
lyte. This leads to a dramatic change of the solution phy- 
sicochemical properties, such as its electrical conductivi- 
ty, at the solute concentration «o,max( T)- 

The estimations performed of the lower N cmin (T,n 0 ) 
and upper N cma JiT,n 0 ) limits of the nucleation critical 
number size N c (T,n 0 ) allow determination of simple 
analytical expressions for the surface tension a(T) and 
the equilibrium bulk energy T) per solute molecule: 


TABLE II. In this table all calculated quantities are obtained at T=298 K. 


No. 

Solute 

a 

J3 

K,mw (rtm+0.1) 

■^c.min^spin ) 

H c , m Jn M +0A) 

N c , max spin ^ 

1 

NaBr 

7.0206 

0.2198 

4.9374 

1 

14.4593 

1 

2 

KCl 

7.1265 

0.1007 

9.4578 

1 

19.8006 

1 

3 

KBr 

6.6947 

0.1466 

6.9080 

1 

16.5334 

1 

4 

NH 2 C1 

6.4823 

0.1662 

6.2668 

1 

16.4945 

1 

5 

(NH 4 ) 2 S0 4 

6.8338 

0.1292 

8.1944 

1 

33.8936 

1 
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0 20 40 60 n 

solute concentration ( in molal units ), n g 


FIG. 5. Dependence of the lower 
^c,min (T,n 0 ) and upper N c>mhX ( T,n 0 ) limits of 
the critical number size N c {T,n 0 ) on the mean 
solute concentration n 0 at 7=298 K for su- 
persaturated aqueous solution of ammonium 
sulfate. 


a(T)=-ln[^ sat (T)] , -ln[£(T)n m (T)] , 

where all concentrations are given in molal units. The 
calculated values for a( T) and T) at 7 = 298 K for the 
investigated supersaturated solutions can be found in 
Table II. 

The conclusions presented above summarize that com- 
bination of the high precision SVELT experiments with 
the theory developed for supersaturated electrolyte solu- 
tions has allowed us to obtain a number of interesting re- 


sults in the physics of supersaturated electrolyte solu- 
tions. 
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